Formation of a narrow baryon resonance with positive 
strangeness in collisions with Xe nuclei 
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' The data on the charge-exchange reaction i^"'"Xe K^pX.e\ obtained with 

' the bubble chamber DIANA, are reanalyzed using increased statistics and updated 

H ' selections. Our previous evidence for formation of a narrow pK^ resonance with 

mass near 1538 MeV is confirmed and reinforced. The statistical significance of 
the signal reaches some Sex (6(t) when estimated as S/^/B {S/ VB + S). The mass 
and intrinsic width of the baryon are measured as m = 1538 it 2 MeV and 
r = 0.39 ±0.10 MeV. 
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The baryons built of four quarks and an antiquark as the lowest Fock component, 
referred to as pentaquarks, are not forbidden by theory and have been discussed for 
many years [Ij. Some properties of light pentaquark baryons forming the lowest SU{3) 
multiplet, the antidecuplet, have been predicted by Diakonov, Petrov, and Polyakov in 
the framework of the chiral quark-soliton model [2]. In particular, they predicted m ~ 
1530 MeV and F < 15 MeV for the mass and width of the explicitly exotic baryon 
with positive strangeness, the Q^{uudds) that should decay to the nK^ and pK^ final 
states. More recent theoretical analyses suggest that the intrinsic width may be as 
small as ~ 1 MeV or even less [3]. Narrow peaks near 1540 MeV in the nK^ and pK^ 
mass spectra were initially detected in low-energy photoproduction [1] and in the charge- 
exchange reaction K^n pK^ |5]. Subsequently, both experiments were able to confirm 
their initial observations ^ [7] . Moreover, increasing the statistics of the charge-exchange 
reaction allowed DIANA to directly estimate the intrinsic width: F = 0.36 ±0.11 MeV 
[7]. Other searches for the 6^ baryon in different reactions and experimental conditions 
yielded both positive and negative results, see the review papers [8] and [9] and references 
therein. The bulk of null results can be probably explained by the extreme smallness of 
the G"*" width that implies the smallness of production cross-sections [TO] . 

The charge-exchange reaction K^n K^p on a bound neutron, that is investi- 
gated by DIANA and BELLE [TTj, is particularly interesting because it allows to probe 
the G"*" intrinsic width in a model-independent manner. The existing data on low-energy 
K^D scattering have been found to leave room for a pK^ resonance with mass near 1540 
MeV, provided that its width is less than 1 MeV [H IH EllISl [IS]. An important 
advantage of the reaction K^n K^p is that the strangeness of the final-state pK^ 
system is a priori known to be positive. In this paper, the DIANA data on the charge- 
exchange reaction K^'K.e —>■ K^pXe' are reanalyzed using increased statistics and updated 
selections. 

The detector and the experimental procedure have been detailed in [7] and refer- 
ences therein. Briefly, the non-magnetic bubble chamber DIANA was filled with liquid 
Xenon and exposed to a separated beam with p = 850 MeV from the 10-GeV proton 
synchrotron at ITEP, Moscow. A traveling through the chamber volume is con- 
tinuously decelerated by ionization, and its momentum at collision point is a tabulated 
function of its path in liquid Xenon. Thereby, we are able to scan K~^X.e collisions over 
an extended interval of momentum between zero and some 700 MeV. On the other 
hand, soft charged secondaries are efficiently identified by bubble sensity and momentum- 
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analyzed by range in Xenon. (Note that this simple and elegant scheme will be imitated 
in future high-statistics studies of the charge-exchange reaction, see llTj-) The estimate 
of momentum based on measured position of the interaction vertex has been verified 
by detecting and reconstructing the tt+tt+tt^ decays in flight. 

Scanning of the film yielded nearly 55 000 events with visible Kg decays, Kg — >■ 



TT ' TT 



and K'^s 



0^0 



TT TC 



in the full fiducial volume of the bubble chamber. These Ko 



decays could be associated with primary fC+Xe vertices with various multiplicities of 
secondary particles. Finally, events with a single proton and K^ tt+tt^ in the final 
state were selected as candidates for the charge-exchange reaction K~^n — > K^p on a 
bound neutron. The selected events are then fully measured and reconstructed in space. 
These measurements are still in progress. The quality of the data is best reflected by 
experimental resolution on effective mass of the pK^ system, estimated as cxm — 3.5 MeV 
by error propagation for observed events. 
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Figure 1: The incident K~^ momentum for measured events of the reaction K~^X.e — 
K^pXe'. The shaded histogram results from restricting the Kg proper lifetime: r < 3tq. 



At present, measurements have been restricted to the region Lx > 520 mm, where 
Lk is the length of the K^ path in liquid Xenon before the collision. (Note that there 
is no one-to-one correspondence between Lk and K~^ momentum, because the original 
beam momentum varied by some ±20 MeV in different exposures.) Shown in Fig. [1] 
is the lab momentum of the incident K'^ for all measured events with Kg and proton 
momenta above 160 and 170 MeV, respectively (instrumental thresholds). Note that the 
statistics of the charge-exchange reaction has been increased by some 11% as compared 
to [7], and that newly measured events are mostly in the high-momentum region. For 
further rejection of Kg mesons that may have scattered by small angles in liquid Xenon 
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but passed the pointback criteria, we then apply the selection r < 3ro where r is the 
measured proper lifetime and tq is its tabulated mean value. 

Figure [2] shows the distributions of measured pKg events with r < Stq in the 
following kinematic variables : 

• laboratory momenta of the Kg and proton ; 

• in the pKg rest frame, the Kg and K'^ angles with respect to the pKg boost from 
the laboratory frame, denoted as Q*{Kg) and Q*{K'^) ; 

• in the laboratory frame, the longitudinal momentum of the pKg system as a whole, 
pdpKl) ; 

• the effective target mass mlf^^, computed assuming binary kinematics. 

Also shown are corresponding distributions of simulated events prior to any proton or Kg 
rescattering in nuclear matterj^ that have been normalized to the number of measured 
events. The Kg and proton momenta are seen to be reduced by rescattering, see Figs. [2^ 
and[2)D. The effect is stronger for protons in agreement with the assumption that the free 
path in nuclear matter is bigger for kaons than for protons [19j. Likewise, the Kg tends 
to travel in the forward hemisphere in the pKg frame, see Fig. [2]3, because rescattering in 
nuclear matter is more probable for protons than for kaons. The K~^ travels exclusively 
in the forward direction cos0*(-ft'^) > for simulated events, see Fig. [21i, but also in 
the backward hemisphere for live events. For the pKg longitudinal momentum plotted 
in Fig. [2^, the simulation predicts that unrescattered events should populate the region 
Pl > 0, whereas the live events extend down to some -300 MeV. The distribution of the 
effective target mass m^f^g much broader for the live data than for simulated events, 

^The non-resonant charge-exchange reaction K^n K'^p in nuclear environment is simulated using 
a simple Monte-Carlo procedure. The total cross-section as a function of collision energy is parametrized 
using the existing data [18] . and the angular distribution in the pKg rest frame is assumed to be isotropic. 
The total energy of a bound neutron is parametrized as E„ = rn„ — e where m„ is the mass of a free 
neutron and e ~ 22 MeV is the mean binding energy for nucleons in the Xenon nucleus. For the same 
nucleus, we use a realistic form of the Fermi-momentum distribution with maximum near 160 MeV/c. 
Pauli blocking for secondary protons is approximated by the cut pp > 170 MeV on proton momentum. 
The flux of incident mesons as a function of momentum at collision point is inferred from 
the observed distribution of range in Xenon before interaction or decay, see [Q. The experimental 
uncertainties and measurements errors are included in the simulation. Rescattering in the nucleus is not 
accounted for. 
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Figure 2: Laboratory momenta of the (a) and proton (b), the cosines of the angle 
(c) and of the angle (d) in the pK^ rest frame, the longitudinal momentum of the 
pK^ system in the lab frame (e), and the effective target mass m!^-^^ (f). The dotted 
histograms are the corresponding distributions of simulated events prior to any proton or 
rescattering in nuclear matter, that have been normalized to the number of measured 
events. 



see Fig. [Sf, because rescattering disrupts the two-body kinematics of the reaction K^n —>■ 

Compared in Fig. [3] are the distributions of two interrelated variables: the effective 
mass of the detected pK^ system, m{pKg), and the center-of-mass energy of the K~^n 
collision assuming a stationary free neutron as the target, ^/s. That the m{pK^) spectrum 
is broader than the ^/s distribution is due to Fermi motion of the target neutron, as well 
as to rescattering of the secondary proton and in nuclear medium that distorts the 
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mass of the originally formed pK^ system. 




pK° mass, GeV ^ assuming free neutron, GeV 

Figure 3: The effective mass of the detected pKg system (a) and the center-of-mass energy 
of the K~^n collision assuming a stationary free neutron as the target (b). The dotted 
histograms depict the distributions of simulated unrescattered events normalized to the 
number of measured events. 

Qualitatively, collisions with bound and free neutrons should form pKg systems 
with similar effective masses (prior to proton and Kg rescattering) . The mass difference 
Am = m{pKg) — ^/s and the asymmetry = {m{pKg) — y/l)/{m{pK^) + y/s) are plotted 
in Fig. m for live and simulated events of the charge-exchange reaction. For simulated 
events, the "original" value of the pK^ mass prior to any rescattering is substituted. The 
observed Am and Am distributions are seen to extend to negative values far in excess 
of the smearing expected from neutron binding in the Xe nucleus. In order to suppress 
rescattered events, we use the selection \Am\ < 0.015 suggested by the simulation. The 
effect of this selection on the distribution of cos Q*{K^), the cosine of the emission 
angle in the pK^g rest frame, is shown in Fig.^t'- By suppressing rescatterings, the selection 
\Am\ < 0.015 is seen to render the latter distribution more uniform. On the other hand, 
this selection brings the m^^j.g distribution to better agreement with the mass of a bound 
neutron, see Fig. [5]d. 

The m{pKg) distribution under the selection \Am\ < 0.015 shows an enhancement 
near 1538 MeV whose width is consistent with instrumental smearing alone, see Fig. [61 
In the fit of this m{pKg) spectrum to a Gaussian on top of a third-order polynomial, the 
width of the Gaussian is constrained to a = 3.5 MeV as estimated for a zero- width reso- 
nance by propagating measurement errors for observed pKg pairs. This fit returns nearly 
100 events above background. In order to further reduce the background from rescatter- 
ings, an additinal selection is used: 6;^ < 100° and 0p < 100° for polar angles of the Kg 
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Figure 4: The mass difference Am = m — \fs (a) and the asymmetry Am = {m — y/s) / (m+ 
^/s) (b) for observed events. Corresponding distributions of simulated events are shown 
in (c) and (d). Shaded areas in (b) and (d) are for \Am\ < 0.015. 




Figure 5: The cosine of the emission angle in the pK^ frame (a) and the effective 
target mass assuming binary kinematics (b). Shaded histograms are for the selection 
\Am\ < 0.015. Dark histograms result from an additional selection Pl{pK^s) ^ ^0 MeV. 



and proton as in [7], or pL{pKg) > 80 MeV as suggested by Fig. [2^, or cos Q*{K~^) > 0.2 
as suggested by Fig. [2]i. Each additional selection further flattens the cos 9cm distribution 
and suppresses the downward tail of the ml^j.^ distribution (for p^ > 80 MeV, this is illus- 
trated in Fig. [5]) . On the other hand, each selection further emphasizes the pK^ peak near 
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Figure 6: Effective mass of the pKg system upon applying tlie selection \Am\ < 0.015 
(a). An additional selection is 6^ < 100° and Qp < 100° in (b), Pl{pK^s) > §0 MeV 
in (c), and cosQ*{K~^) > 0.2 in (d). Depicted by dotted curves are fits to a Gaussian 
on top of a third-order polynomial. The width of the Gaussian has been constrained 
to o" = 3.5 MeV. The solid dots are simulated distributions assuming no rescatterings, 
each normalized to the corresponding experimental distribution by area outside the peak 
region 1529 < m{j)K^) < 1547 MeV. The simulated mass spectra have been scaled by a 
factor of 0.5 in order to avoid confusion with the fitting function. 



1538 MeV, see Figs. [6]d, [Ht, and[6li. The statistical significance of the signal reaches some 



8a {6a) when estimated as S/yB {S/y/B + S). Also shown in Fig. [6] are the simulated 
distributions of the pKg effective mass prior to any rescatterings under similar selections. 
Each of these has been normalized to the corresponding experimental distribution by area 
outside the peak region 1529 < m{pKg) < 1547 MeV. That peak position is very close to 
the maximum of the non-resonant background may lead to overestimating the signal. So 
we additionally fit the difference between the experimental and simulated distributions, 
see Fig. [71 The fits of subtracted mass spectra return virtually the same signals as in 
Fig. [HI suggesting that the fitting procedure is robust. 

In ref. [7] where the selections 9^ < 100° and 0p < 100° were used, we argued 
that the bulk of the G"*" signal should arise from incident mesons with momenta in 
a restricted interval. Plotted in Fig. [Hj is the pK^ effective mass under the selections 
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Figure 7: The difference between the experimental and simulated spectra of the pK\ 
effective mass, as shown in Fig. [61 



as in Fig. [6b plus the additional selection 445 < p{K^) < 525 MeV as in [7]. This 
is separately done for the complete sample of the charge-exchange reaction and for the 
part of the sample analyzed in [7]. The signal in Fig. [8)d is virtually the same as that 
reported in [7]. Combining the selections \A„i\ < 0.015 and 445 < p{K^) < 525 MeV is 
of course a redundant procedure used only for comparing with our previous results. That 
the selection \Ajn\ < 0.015 is not tuned to any particular value of the pK^ effective mass 
is its obvious advantage. 

Intrinsic width of a resonance formed in an s-channel reaction like K^n K^p 
can be estimated by comparing the signal magnitude with the level of non-resonant back- 
ground under the peak, see e.g. [13J. However, this method cannot be directly applied to 
collisions with heavy nuclei because of secondary interactions of the and proton 
in nuclear medium. That is, the non-resonant background under the peak is an unknown 
mixture of unrescattered and rescattered events, whereas a true signal should consist of 
unrescattered events only. As soon as the B"*" decay width is ~ 1 MeV or less, the peak 
will not be depleted by the and proton rescatterings because the bulk of produced 9+ 
baryons will decay upon leaving the nucleus. The corresponding "rescattering-free" back- 
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Figure 8: The pKg effective mass under the selections as in Fig. [Ht" (|^m| < 0.015, 
Qk < 100°, and Qp < 100°) plus the additional selection 445 < p{K+) < 525 MeV. 
The data for the complete sample of the charge-exchange reaction and for the part of the 
sample analyzed in [7] are shown in (a) and (b), respectively. 



ground can only be estimated by simulating the effective mass of the originally produced 
Kg and proton prior to any rescatterings, rriolpKg). 

For the simulated charge-exchange collisions K^n K^p on a bound neutron in 
the bubble chamber DIANA, the momentum at collision point and the "undistorted" 
mass mo{pKg) are plotted in Fig. [9l The Kg lab momentum is restricted to px > 140 
MeV which is the (effective) threshold for detecting Kg —>■ tt+tt^ decays in the scan. 
The upper histograms are for all collisions in the full fiducial volume, i.e., over the full 
range of Lk — the K~^ track length in Xenon before the collision. The shaded area of the 
K~^ momentum distribution corresponds to the restricted fiducial volume > 520 mm 
where throughput measurements were done. The absolute normalization of the simulated 
distributions of Fig. [9] is based on the scanning information. In the region > 520 mm, 
the scan found 7900 ±500 events with a Kg tx^tx~ candidate in the final state. Of these 
events, (60 ± 7)% are estimated to survive upon rejecting K^ mesons with L < 2.5 mm, 
protons and Kg mesons that have reinteracted in liquid Xenon, and unmeasurable events. 
The resultant number of events, 4740 ±630, is then used for absolute normalization of the 
simulated i^^-momentum distribution for > 520 mm, see the shaded area in Fig. [9^. 
The distribution of the same events in mo{pKg) is shown by the corresponding shaded 
histogram in Fig. [9)3. And finally, the absolute level of the non-resonant background can 
be estimated upon applying the aforementioned selections to simulated events: px > 160 
MeV, Pp > 170 MeV, r < 3tq, and 1^4^,1 < 0.015. The resultant niQ^pKg) distribution, 
depicted by the dark-shaded histogram in Fig. [Hb, features 705 ± 94 events in the mass 
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Figure 9: The momentum (a) and the "undistorted" mass mQ{pK'g) (b) for simulated 
charge-exchange colhsions K^n K^p in the bubble chamber DIANA. The upper his- 
tograms are for all collisions with px > 140 MeV in the full fiducial volume. The shaded 
histograms are for track lengths Lk > 520 mm (the region of throughput measure- 
ments). The dark-shaded histogram in (b) results from the adopted selections pk > 160 
MeV, Pp > 170 MeV, r < Stq, and \Am\ < 0.015. The open-white corridor in the latter 
histogram depicts the mass region 1529 < rriolpKg) < 1547 MeV. The normalization is 
explained in the text. 

region 1529 < mo(pA'g) < 1547 MeV. 

Using the observed signal and simulated background as input and assuming J = 
1/2, we then estimate the intrinsic width of the pK^ resonance as [13] 

^ _ NP^^'^ ^ a^^ ^ Am 
- ATbkgd ^ i07mb ^ BiBf' 

where N'^'^^^ and A^'^'^^d ^^-g ^]-^g numbers of events in the peak and in the background 
under the peak, cr*"^ = 4.1 ± 0.3 mb is the measured cross section of the charge-exchange 
reaction K'^n -—>■ K^p [T8l|, Bi and Bj are the branching fractions for the initial and final 
states {Bi = Bf = 1/2 for either 1 = and / = 1), and Am is the m{pK^) interval 
under the peak that is populated by N^^^'^ background events. Selecting a mass interval 
1529 < m{pK^) < 1547 MeV and substituting A^p^^"^ = 99.9 ±18.2 events from Fig. Eli and 
^bkgd ^ 705±94 events from Fig. [9b, we obtain F = 0.39±0.10 MeV where the error does 
not include the systematic uncertainties of the simulation procedure. A systematic shift of 
up to +20% may result from taking into account the sequential reaction K^N K^N, 
K^n K^p [20]. The value of F derived in this analysis is consistent with our earlier 
estimate [7] and with the upper limit set by BELLE [1 1] . 

To conclude, we have reanalyzed the DIANA data on the charge-exchange reaction 
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K~^X.e — i> K^pXe' using increased statistics and updated selections. We have obtained 
new strong evidence for formation of a pentaquark baryon with positive strangeness in the 
charge-exchange reaction K~^n —>■ K^p on a bound neutron. The statistical significance 
of the signal reaches some 80" (6cr) when estimated as S/\fB {S/y/B + S). The mass and 
intrinsic width of the G"*" baryon are measured as m = 1538 ± 2 MeV and F = 0.39 ± 0.10 
MeV. 

We wish to thank Ya. I. Azimov, L. N. Bogdanova, and I. I. Strakovsky for useful 
comments. This work is supported by the Russian Foundation for Basic Research (grant 
07-02-00684). 
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